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In this paper we use a traveling wave reduction or a so-called spatial approximation to comprehen- 
5_( ■ sively investigate the periodic solutions of the complex cubic-quintic Ginzburg-Landau equation. 

CLT The primary tools used here are Hopf bifurcation theory and perturbation theory. Explicit results 

, are obtained for the post-bifurcation periodic orbits and their stability. Generalized and degener- 

ate Hopf bifurcations are also briefly considered to track the emergence of global structure such as 
homoclinic orbits. 
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I. INTRODUCTION 



The cubic complex Ginzburg-Landau equation (CGLE) is the canonical equation governing the weakly 
nonlinear behavior of dissipative systems in a wide variety of disciplines (Tlj . In fluid mechanics, it is also 
often referred to as the Newell- Whitehead equation after the authors who derived it in the context of Benard 
convection [Til [l6| . 

As such, it is also one of the most widel y s tudied nonlinear equations. Many basic properties of the 
equation and its solutions are reviewed in (4l.llfjj. together with applications to a vast variety of phenomena 
including nonlinear waves, second-order phase transitions, superconductivity, superfluidity, Bose-Einstein 
condensation, liquid crystals and string theory. The numerical studies by Brusch et al. [e| [9| which primarily 
consider periodic traveling wave solutions of the cubic CGLE, together with secondary pitchfork bifurcations 
and period doubling cascades into disordered turbulent regimes, also give comprehensive summaries of other 
work on this system. Early numerical studies [U [22[ and theoretical investigations [lH [26[ of periodic 
(f") ■ solutions and secondary bifurcations are also of general interest for our work here. 

Certain situations or phenomena, such as where the cubic nonlinear term is close to zero, may require the 
inclusion of higher-order nonlinearities leading to the so-called cubic-quintic CGLE. This has proved to be a 
rich system with very diverse solution behaviors. In particular, a relatively early and influential review by van 
Saarloos and Hohenberg [33( , also recently extended to two coupled cubic CGL equations @, |32[ , considered 
phase-plane counting arguments for traveling wave coherent structures, some analytic and perturbative 
solutions, limited comparisons to numerics, and so-called "linear marginal stability analysis" to select the 
phase speed of the traveling waves. 

Among the multitude of other papers, we shall only refer to two sets of studies which will directly pertain 
to the work in this thesis. The first class of papers 0, HH, H3] and [H,[l7| used dynamical systems techniques 
to prove that the cubic-quintic CGLE admits periodic and quasi-periodic traveling wave solutions. 

The second class of papers @, primarily involving numerical simulations of the full cubic-quintic CGL 
PDE in the context of Nonlinear Optics, revealed various branches of plane wave solutions which are referred 
to as continuous wave (CW) solutions in the Optics literature. More importantly, these latter studies also 
found various spatially confined coherent structures of the PDE, with envelopes which exhibit complicated 
temporal dynamics. In [2j, these various structures are categorized as plain pulses (or regular stationary 
solutions), pulsating solitary waves, creeping solitons, slugs or snakes, erupting solitons, and chaotic solitons 
depending on the temporal behavior of the envelopes. In addition, note that the speed of the new classes 
of solutions may be zero, constant, or periodic (since it is an eigenvalue, the speed may be in principle also 
quasiperiodic or chaotic, although no such cases appear to have been reported). All indications are that 
these classes of solutions, all of which have amplitudes which vary in time, do not exist as stable structures in 
Hamiltonian systems. Even if excited initially, amplitude modulated solitary waves restructure into regular 
stationary solutions @ . Exceptions to this rule are the integrable models where the pulsating structures are 



nonlinear superpositions or fundamental solutions [30|. Hence, these classes of solutions are novel and they 
exist only in the presence of dissipation in the simulations of [2]. Also, secondary complete period doubling 
cascades of the pulsating solitons leading as usual to regimes of chaos are also found. This last feature for 
numerical solutions of the full cubic-quintic PDE is strongly reminiscent of the period doubling cascades 
found in [1, Q for period solutions of the traveling wave reduced ODEs of the cubic CGLE. 

Motivated by the above, we begin a fresh look at the traveling wave solutions of the cubic-quintic CGLE 
in this paper. Besides attempting to understand the complex numerical coherent structures in [2J, one other 
goal is to build a bridge between the dynamical systems approach in [ID, [H], H^-flll G3 and the numerical 
one in 0,|3l|. Given the importance of the cubic-quintic CGLE as a canonical pattern-forming system, this 
is clearly important in and of itself. However, a word of warning is in order here. Some of the features in 
may well be inherently spatio-temporal in nature, so that a spatial traveling-wave reduction may not be 
sufficient to completely capture all aspects. Indeed, there is some evidence along these lines [§[. 

In this paper we begin by using a traveling wave reduction or a so-called spatial approximation to investi- 
gate the periodic solutions of the CGLE. The primary tool used here is Hopf bifurcation theory. Immediately 
following the Hopf bifurcations we construct the periodic orbits by using the method of multiple scales. 

The remainder of the paper is organized as follows. We first analyze the stability of fixed points in Section 
2 and the onset of instability via a Hopf bifurcation, which may be either supercritical or subcritical. Then 
stability of periodic orbits is presented in Section 3 where we derive analytical expressions for the periodic 
orbits resulting from this Hopf bifurcation, and for their stability coefficients, by employing the multiple scales 
method. Section 4 considers numerical solutions and discusses the results. Generalized and degenerate Hopf 
bifurcations have also been considered to track the emergence of global structure such as homoclinic orbits 
in Section 5. 



II. STABILITY ANALYSIS OF FIXED POINTS 



In this section, we conduct a stability analysis of individual plane wave solutions using regular phase plane 
techniques. This was already done for the alternative formulation of the traveling wave ODEs given in [33j . 
We provide a brief derivation for our system (|II.4a[) . (|II.4b[) and pi.4cl) for completeness and future use. 
However, a much more complex question is the issue of categorizing and elucidating the possible existence 
of, and transitions among, multiple plane wave states which may co-exist for the same parameter values 
in pi.ip (corresponding to the same operating conditions of the underlying system ). S uch behavior is well- 
documented in systems such as the Continuous Stirred Tank Reactor System [l|, |30| . For a system such 
as pi.l[> and the associated ODEs pi.4a[) . pi.4b|) and pi.4cl) . the large number of parameters makes a 
comprehensive parametric study of co-existing states bewilderingly complex, if not actually impracticable. 

We shall consider the cubic-quintic CGLE in the form [33| 

d t A = eA + (bi + icx)dlA - (b 3 - ic 3 )\A\ 2 A - (b 5 - ic b )\A\ A A (II.l) 

noting that any three of the coefficients (no two of which are in the same term) may be set to unity by 
appropriate scalings of time, space and A. 

For the most part, we shall employ the polar form used in earlier treatments 0, [§3[ of the traveling wave 
solutions of (III.ll) . This takes the form of the ansatz 

A(x,t) =e- iut A(x-vt) 

= e- lut a{z)e l ' t ' {z) {11.2) 

where 



z = x — vt (II. 3) 

is the traveling wave variable and ui and v are the frequency and translation speed (and are eigenvalues). 



Substitution of pL2|) / (|TT73 ]) in pTT]) leads, after some simplification, to the three mode dynamical system 

a z = b (II. 4a) 

b z = aijj 2 — 71 72a + v(bib + ci^paj — 73a 3 — 74a 5 

7i 75 + v {~ ~ ~ 76a 2 - 77a 4 



2# 

Wz = 

a 



(II.4b) 
(II.4c) 



where ip = <f> z . Note that we have put the equations into a form closer to that in [9], rather than that in 
[33l | , so that (|II.4|1 is a generalization of the traveling wave ODEs in [H to include the quintic terms and the 
constant terms 71 — 77 are given as functions of the system parameters in the following manner: 
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— b\LU + c\e 
bic 3 + C1&3 
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From (|H.2[) . a fixed point (00,0,^0) 01 (|H-4[) corresponds to a plane wave solution 

A(x,t) = aQe i{^oz- u t)+ie 



(II.5) 



with 6 an arbitrary constant. 

The fixed points of pi.4a[) . (|II.4b[> and (|II.4c|) may be obtained by setting b = (from ||II.4a|) ) in the right 
hand sides of the last two equations, solving the last one for tjj, and substituting this in the second yielding 
the quartic equation 



0L2X + ot\x + eta = 



(II, 



with 
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(II. 7a) 
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= 71(73+ h 


> 2 7 5 7 6 
' b\v 2 


(II.7e) 
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75 71/, . \ 
= h2 2 ^ (0l72 +C175J • 


(II.7f) 



Thus, with ao = \fx for x any of the four roots of (|II.6I) . we have a plane wave solution of the form (|II.5|) . 
For each of the four roots Xi , i — 1 , . . . , 4 of (|II.6|) corresponding to a fixed point of (|II.4a[) , (|II.4bj) and piTcft 
or a plane wave yjxl e'W 2_,l ' i )+ ,Ji 1 the stability may be determined using regular phase-plane analysis. The 



characteristic polynomial of the Jacobian matrix of a fixed point Xi 
be expressed as 



of (|II.4aj) . pi.4bl) and piTcft may 



A 3 + <5iA 2 + 5 2 A + S3 = 



(II. 



where 



4 



5j = 2M 7l (II.9a) 

5 2 = 3^ 2 + 71 [72 - a 2 (3 73 + 5a 2 74 ) - v(3c^ - v)] (II.9b) 

5 3 = -2a 2 7i(7 6 + 2a 2 j 7 )(-2ip + Ci 7 iu) 

+ &i7i«[-'0 2 + 7i(72 - 3a 2 7 3 - 5a 4 7 4 + Citpv)} (II.9c) 



where the fixed point values (ai,ipi) = {^fxl^Ji) are to be substituted in terms of the system parameters. 
Note that ipi is obtained by setting a = a,i = y/~Xi, and b = in the right side of (III.4cp . 

For (ao,0,?/>o) to be a stable fixed point within the linearized analysis, all the eigenvalues must have 
negative real parts. Using the Routh-Hurwitz criterion, the necessary and sufficient conditions for (|II.8|) to 
have i?e(Ai ; 2,3) < are: 

81 > 0, 5 3 > 0, 5^2 - S 3 > 0. (11.10) 

Equation (jII.10[) is thus the condition for stability of the plane wave corresponding to Xi . 

On the contrary, one may have the onset of instability of the plane wave solution occurring in one of two 
ways. In the first, one root of (|II.8|) (or one eigenvalue of the Jacobian) becomes non-hyperbolic by going 
through zero for 

S 3 = 0. (11.11) 

Equation (|II.11[) is thus the condition for the onset of "static" instability of the plane wave. Whether this 
bifurcation is a pitchfork or transcritical one, and its subcritical or supercritical nature, may be readily 
determined by deriving an appropriate canonical system in the vicinity of (|II.11|) using any of a variety of 
normal form or perturbation methods 0, EH EH . 

One may also have the onset of dynamic instability ( "flutter" in the language of Applied Mechanics) when 
a pair of eigenvalues of the Jacobian become purely imaginary. The consequent Hopf bifurcation at 

SiS 2 -6 3 = (11.12) 

leads to the onset of periodic solutions of pi.4al) , pi.4b[) and (III.4c[) (dynamic instability or "flutter" ) . These 
periodic solutions for a{z) and ip{z), which may be stable or unstable depending on the super- or subcritical 
nature of the bifurcation, correspond via (|II.2j) to solutions 

A(x, t) = a{z)e t( 5 (11.13) 

of the CGLE (jll.lj) which are, in general, quasiperiodic wavetrain solutions. This is because the period of ip 
and ui are typically incommensurate. Eq. (|II.13|) is periodic if w = 0. 



III. STABILITY ANALYSIS OF PERIODIC ORBITS 



In this section we will use the method of multiple scales to construct analytical approximations for the 
periodic orbits arising through Hopf bifurcation of the fixed point of the CGLE equation. For the systems of 
differential equations given by pi.4a[) . (jll.4b|) and (|II.4cj) . the physically relevant point is given by (ao, 0, t/>o) 
where ipi is obtained by setting a = cii = ifxl, in 

t/H = ^-^(76 + a 2 77) (nL1) 



and Xi is one of the roots of the fixed point equation (|II.6|) . We will choose the parameter e which represents 
the linear gain or loss as the control parameter. The limit cycle is determined by expanding about the fixed 



point using progressively slower time scales. The expansion takes the form 

3 



a = a + J2 ^ a n(Zo, Z lt Z 2 ) + ■ ■ ■ , (III.2) 

n=l 
3 

b = B Q + Y / S n B n {Z Ql Zi, Z 2 ) + ■ ■ ■ , (III.3) 

71=1 

3 

^ = ^o + Y, Sn MZo, Z lt Z 2 ) + >>-, (III.4) 



where Z n = 5" t and 5 is a small positive non-dimensional parameter that is introduced as a bookkeeping 
device and will be set to unity in the final analysis. Using the chain rule, the time derivative becomes 

^ = D + SDi + 5 2 D 2 + ■■■ , (HI.5) 

where D n — d/dZ n . The delay parameter e is ordered as 

e = e + <5 2 e 2 , (111.6) 

where eo is the critical value such that (III.10[) is not satisfied, (i.e. eo is a solution of pi.l2[0 . This is standard 
for this method, as it allows the influence from the nonlinear terms and the control parameter to occur at 
the same order. 

Using (|III.2[) - (IIII.6I) in (|II.4aj) — (|II.4cj) and equating like powers of 5 yields equations at 0(6 t ), i = 1,2, 3 
of the form: 

Li(ai,Bi,tpi) = S^i, (III.7) 
L 2 (a i ,B i ,i> i ) = S it2 , (III.8) 
L a (a i ,B i ,^ i ) = S i<3 , (III.9) 

where, the Li, i = 1,2,3 are the differential operators 

L 1 (a i ,B i ,i/j i ) = D ai-B i = Si,x, (111.10) 
L 2 (a i ,B i ,ip i ) = D Bi - i/j 2 ai - 2a ip ^i 

+ 7i{72oai + v[b 1 B i +ci(^ a,: + a ipi)] 

- 3"/3oalcii - 574oaoai} = 5i,2, (III. 11) 
L 3 (ai,Bi,ij)i) = ao(A)V'j) + 2(ip Bi + B ipi) 

- 7i{75oa* + v[c 1 B l - bi(ifj a,i + a ipi)] 

- 376oao a i _ 5770agai} = S i>3) (III. 12) 

where -f p = 7 p0 + S 2 7 p2 with p = 2, 3, • • • ,7, the source terms Sjj for i, j = 1, 2, 3 at 0(5), 0(<5 2 ), and 0(6 3 ) 
are given by the following: 

0(6): 



Si,i = o 


(III. 13) 




(111.14) 


#1,3 - 0. 


(111.15) 



0(6 2 ) : 



S 2A = -Dioi (111.1$ 
£2,2 = -DxBx + a ^l + 2ip aiipi - 71(72200 + vciipiai - 3~f 30 a al 

- 732«o - !074oaoai - 742^) (111.17) 
5 2: 3 = -a £>i^i - ai-DoVh - 2ipiBi + 71 [(75200 - v&i^iai) 

- (376oaoa? + 762ao) - (lC^oa^ + 772flo)]- (III. 18) 

0(6 3 ) : 

S 3 ,i = -D ia2 - D 2 a x (111.19) 
•53,2 = —D\B 2 — D 2 Bi + 2a 2 ip ipi + 01(2-00^2 + i>i) + 2ao4>iip2 

- 71(72201 + vci(ipia 2 + ip2di) - [73o(a? + 6a aia 2 ) + 37320^1] 

- [ 74o (10aga? + 20a 3 a ia2 ) + 57 42 ao a i]} (111.20) 
£3,3 = -D^j-2 - D 2 ipi - a^Dxipi + D ip 2 ) - a 2 D ip x - 2(^ 1 B 2 + ip 2 B x ) 

+ 7i {752ai - vbi(^!a 2 + ip 2 ai) - [76o(a? + 6a aia 2 ) 

+ 3762flo a i] - [77o(10aoa? + 20a 3 ,aia 2 ) + 5772^1]}. (111.21) 
Also, (|III.7p may be solved for Bi in terms of a^. and ipi. Using this in (|III.8|) yields ipi 

0i = (HI.22) 

where 

Oi = -DoSi,i + D a,i - tp + 7i{72oa-i - 37300^1 - 574oflo a * 

+ + Doch) + ci]} - Si, 2 (111.23) 

and 

<j>\ = 2oo^o - t'7icia . (III. 24) 

Using (|III.22j) and the equation for Bi in (|III.10[) yields the composite equation: 

L c ai = Ti, (111.25) 

where 

^ = S it a - ^(DoCi) - ^Ci - livhao^- + (2^o - 71^1)^,1, (111.26) 

Q = -D S iA - TivbiS it i - S it2 . (IH.27) 

We shall now use (|III.26I) and (|III.27|) to systematically identify and suppress secular terms in the solutions 
of piI.10|) , piI.ll|MIII.12|t . Let us now turn t0 finding the solutions of 1 0|) . (IIII. 1 1 1) . pil. 1 2 1) . In what 
follows , we shall detail the solution of the above system of equations for the case eo = eoi. In order to 
achieve that we must find first the fixed points. The characteristic polynomial of the Jacobian matrix of a 
fixed point of pi.4a[l . (III.4b[) . pT4c| may be expresses as 

A 3 + 5i\ 2 + S 2 X + S 3 = 0, (111.28) 

as in pi.8[) . and the fixed point values (ai,ipi) are to be substituted in terms of the system parameters. 

The condition 6162 — 63 — yields an involved equation in e which actually can be solved easily numerically 
for eo by the root method . 

For 0(6) the Eqns. (|III.13I) (|III.15|) give SVi = Si, 2 = 5^3 = 0, and hence we may pick a solution for the 



first order as 



ai = a(Z u Z 2 )e MZ + /3(Z u Z 2 )e^ z + 7 (Z lj Z 2 )e 



(111.29) 



where /? = a is the complex conjugate of a and A2 = Ai. As evident for the Routh-Hurwitz condition, 
the a and /3 modes correspond to the center manifold where A^a are purely imaginary and where the Hopf 
bifurcation occurs, while 7 corresponds to the attractive direction or the stable manifold. Since we wish to 
construct and analyze the stability of the periodic orbits which lie in the center manifold, we should take 
7 = so (|III.29|) becomes 



01 = a{Z ll Z 2 )e lu)Z + p{Z 1 ,Z 2 )e-^ z . 
Using (|III.13[) 15|) for i — 1 then the first order fields (ai, V>i) are 

B 1 = D oai = iuae iuZ - iu$e~ iMZ \ 

and also pil.221) becomes 

1 



(111.30) 
(111.31) 



- V^o + 7i (720 + vcx^o - 37 30 aQ - 57 40 a^ 



x foe*"* + Pe^ z ) + ^ 



icuae 



iujZ 



-iujZ 



(111.32) 



Now that the first order solutions pil.30[t - (jlll.32[) 

are known, the second order sourcss S211 S221 $23 may 
be evaluated via (IIII.16p - pII.18p . Using these sources in ([III.26P we obtain T2 which may be written as 



r 2 = r 2 



(0) 



ifW +lfV iwZ + c.c. 



(111.33) 



Setting the coefficients of the secular e l terms (which are the solutions of the homogeneous equation for 



1) to zero, i.e. T 



(i) 



yields 



^ da 



0Z1 



(111.34) 

Using pil.34|) . the second order sources, and assuming a second-order particular solution for a 2 of the form: 

(111.35) 



a 2 — -\- o> 2 e , 



having the standard form of a DC or time-independent term plus second-harmonic terms, the composite 
equations (IIII.25p - piI.27l) for i = 2, yield 



L c d2 — 1^2 



(0) 



r (2)2iujZ 
1 2 e > 



(111.36) 



which will be solved for the particular solution a 2 °^ , and a 2 by equating both sides of the expression (|III.36[) . 



In terms of the operator L c which is obtained from (|III.26[) , the particular solution takes the form: 

[oo(«ci7i - 2^o )] |2So[ - V'o + 71(720 - 3^730 - 5ag740 + vciipo)] 
ao7i{(750 - 3ao7 60 - 5^770) (vcm - 2ip ) 

vh [71 (720 - 3ao73o - 5ao74o) + V'o] } } - (III. 37) 



4 2) = ~r^ 2) [oo(«ci7i - 2^o)] j6ao-B 7i73o + 10ao- B o7i740 
+ 3ag7i [73o(2zw + vbiji) + 7 60 (vci7i - 2ip )] 



+ 5a^7i [740 {2iuj + vbiji) + 770(^171 - 2ip )] 

+ 2B [Auj 2 - 2ww6i7i + ip 2 - 71(720 + vciipo)] 

+ ao{8iu) 3 — 2iv 2 uib 2 j 2 — 6iuipQ + 27i^q(3wo;ci + 750) 

+ 7i[-2iw(w 2 Ci7i + 720) - WC171750] - w6i7i(-8w 2 +71720 + ^0)}} ■ (III. 38) 

Using (|III.22|) . the second order sources and the equation for Bi in (|III.10[) with i = 2 then we can find 
the second order fields B 2 and -02- Substituting them into the (|III.19[) - (IIII.21[) we find the third order 
sources and we may evaluate the coefficients of the secular term e lut in the composite source 1^ of (|III.26p . 
Suppressing again the secular terms to obtain uniform expansions yields the final equation for the evolution 
of the coefficients in the linear solutions 30[) - (jIII.32[) on the slow second-order time scales 

— = S ia 2 p + S 2 a. (111.39) 
aZ 2 

Writing a — ^Ae te and separating (|III.39I) into real and imaginary parts, yields 

J| = ^W (111.40, 

where Si r and S 2r represent the real parts of Si and S 2 respectively. In the usual way, the fixed points of 
Clianj, Ui. ^2,3) where 

Ai = 0, 



A 2 , 3 = ±2J-^ (111.41) 

give the amplitude of the solution a = ^Ae 10 , with ^2,3 corresponding to the bifurcation periodic orbits. 
Clearly A 2 ^ are real fixed points whenever 

< 0, (IH.42) 

and the Jacobian of the right hand side of (|III.42[) evaluated at A 2 3 is J\a 2 3 — ~S 2r . Clearly, a necessary 
condition for stability is to have S 2r > 0, and for instability S 2r < 0. Thus, the system undergoes: 

a. supercritical Hopf bifurcations when 

S 2r > 0, Si r < 0, (111.43) 



b. subcritical Hopf bifurcations when 

S 2r < 0, Si r > 0. 



(111.44) 











20 



40 



GO 



SO 



100 



z 



FIG. 1: Stable periodic oscillations on the limit cycle 



IV. DISCUSSION OF RESULTS 



In this section, we consider the numerical results which follow from the analysis in the previous section. The 
fixed point equation ()II.6|) can be solved analytically for each fixed point Xi using the program Mathematica, 
for i = 1, • • • , 4. The characteristic polynomial of the Jacobian matrix of a fixed point of (|II.4a|) . (|ll.4b[) . (|II.4cj) 
may be expresses as is (|III.lj) . Since all coefficients on, for i = 1- • • 3 depend on the nine system parameters, 
we fix h = 0.08, 6 3 = -0.65, b 5 = 0.1, c x = 0.5, c 3 = 1, c 5 = -0.07, u = 0, and v = 0.01. The 
possibility of bounded chaotic solitons depends on the system being fairly strongly dissipative near the 
fixed points (ao,0,?/>o) m a significant part of the phase space, with the strong dissipativity ruling out 
the appreciable volume expansion associated with an attractor at infinity, as well as volume-conserving 
quasiperiodic behavior. The trace of the Jacobian matrix for this sets of values at the fixed point (ao, 0, tpo), 
which gives the local logarithmic rate of change of (a,b,ip) phase-space volume V is = J(ao, 0, -0o) = 
—0.0062, so we may anticipate that the orbits may go to an attractor at infinity, since the dissipation is 
weak. 

The four fixed points can be analytically found as a function of only one parameter, in our case we chose 
e as being the free parameter. By choosing "the right fixed point", the Hopf curve a\a 2 — «3 = may be 
solved numerically for e, which gives eo = —0.0000807. The idea is to find the "right" e which will give rise 
to the condition for Hopf bifurcation, (i.e. a± > 0, a 2 > 0, «3 > and ot\a 2 — CV3 < 0). 



We obtain a x = 0.006, a 2 = 0.001, a 3 = 0.0001 and - a 3 = -1.01 10~ 6 for an e = -0.008, i.e. 



Now we will analyze the multiple scales method to construct the analytical approximations for the periodic 
orbits arising through the Hopf bifurcations of the fixed point. The delay parameter e (or the bifurcation 
parameter) is ordered as e = eq + S 2 62, where eo = —0.0000807, and £2 = —0.1. This method allows the 
influence from the nonlinear terms and the control parameter to occur at the same order. For the system 
parameters chosen above, at the fixed points, we get (ao,0, ^0) = (0.0121663,0,-0.00514). 

From (|III.25j) and by the method presented in Section 2, the final equation for the evolution coefficients 
in the linear solutions, on the slow second-order time scale is 



e < e . 



da(Z 1 ,Z 2 ) 
dZ 2 



S 1 a 2 (Zi,Z 2 )(i{Z u Z2) + S 2 a(Zi,Z 2 ), 



(IV.l) 



where Si = -3235.55 + 295.279i and S 2 = 297.074- 32. 26i. Since S 2r = Re(S 2 ) > 0, and S lr = fle(Si) < 0, 
then this situation will correspond to a supercritical Hopf bifurcation. Figure [T] shows the time behavior for 
a(Z) for e = —0.00008 (the supercritical regime). Note that, as anticipated from before, there is a stable 



limit cycle attractor at e, the solution remains positive and bounded while it stays periodic. ^ 
Clearly, similar stable periodic solutions may be obtained for many other parameter sets. For each case, 
the overall solution A(x, t) of the CGLE is, via pi. 131) . a quasiperiodic solution. 

One may also use the above approach to both explain, and extend, the numerical treatment by Brusch 
et al [1, Q of the periodic traveling waves of the CGLE using the bifurcation software AUTO. However, the 
solutions in Brusch et al do not appear to be clearly correlated to the dissipative solitons of the CGLE in 
Akhmediev at al [l[ . Hence, we shall move on next to briefly consider possible generalizations of the above 
treatment. 



V. REMARKS ON GENERALIZED HOPF BIFURCATIONS AND EMERGENCE OF GLOBAL 

STRUCTURE 

One may pursue the line of inquiry based on the traveling waves or spatial ODEs even further to track the 
emergence of global structure. We have done preliminary work along these lines. However, although there 
is a well-established roadmap and it has been implemented in detail for the well known Continuous Stirred 
Tank Reactor System [la . |29| , we are not convinced of its relevance to the actual numerical simulations of 
dissipative solitons P, Hf- Hence, we present it here as a possible future direction to pursue. 

For completeness, let us first consider more degenerate cases where more than one root of the Jacobian is 
non-hyperbolic. In such cases the non-hyperbolic eigenvalues of the Jacobian matrix, may consist of either: 

a. a double zero: Ai,2 =0 A3 € 5i 

b. one zero and a complex conjugate pair: Ai = A2 = A3 

c. a triple zero: Ai,2,3 = 

For the above situations, we have the following sub-cases of the so-called "degenerate Hopf" (HI) bifurcation. 
Each sub-case is given a name: 

Fi: Ai l2 =0,0 (V.l) 
F 2 : Ai, 2)3 = ±iw , (V.2) 

Gi: Ai,2, 3 = 0,0,0 (V.3) 

In these cases, [1 GHH and HHHIH, these (HI) bifurcations may lead to global structure including 
homoclinic orbits, invariant tori, and period doubling to chaos at the (HI) points. One may also work 
perturbatively [f| near these (HI) points as done by Keener for the well-known Continuous Stirred Tank 
Reactor problem. 

Two other degenerate/generalized Hopf bifurcation scenarios are possible. As seen in Chapter 3 pil.401) . 
the normal form for the Hopf bifurcation may be written as 

r = r[a([i) + Cl ([i)r 2 + c 2 (ti)r i + ...] (V.4) 

6 = Lu a + 0(v,r 2 ) (V.5) 

where we have made the identification A — > r, Si r /A — > ci, S% T — > a, and higher order nonlinear terms are 
included. 

The first kind of possible degeneracy (the (H2) kind) occurs if 

a = a' = ... = a {k) = 

a {k+1) ? 0. (V.6) 

This is the so-called k th order (H2) degeneracy and it gives rise to multiple Hopf points and multiple periodic 
orbits. The resulting structure is thus similar to that resulting from a regular Hopf bifurcation, and much 
less complex than the structure produced by (HI) bifurcation. 



A second possible degeneracy in the normal form (|V.4[) corresponds to 
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Ci = c 2 = 

Cm+l 7^ 0. 



(V.7) 



This so-called m th order (H3) degeneracy results in isolated branches of periodic solutions unconnected to 
the main branch. 

When the k th order (H2) degeneracy and the m th order (H3) degeneracy occur simultaneously, the normal 
form (IV.4I) may be rescaled to the form: 



This is the so-called H m k degeneracy. 

In the case of the (H2) degeneracy, the complex conjugate eigenvalues ±iuj at the Hopf point cross the 
imaginary axis tangentially leading, after additional analysis, to multiple periodic orbits. 

For (if 3) degeneracy, one may obtain isolated branches (isolas) of periodic orbits unconnected to the main 
branch. 

However, of greatest interest are the (HI) bifurcations where the Jacobian has more than one non- 
hyperbolic eigenvalue and global structure emerges. These will be pursued in future work, together with the 
(H2) and (H3) cases. 
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